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Abstract 

Two new pulse shapes for communications are presented. The first 
pulse shape is ISI-free and identical with the interpolating function (or 
ISI-free kernel) of a reconstruction formula in shift-invariant spaces 
with Gaussian generator. Several closed form representations in time 
and frequency domain are given including one for an approximation 
that is particularly simple. The second pulse shape is the root of 
the former and obtained by spectral factorization. As a consequence, 
shifted versions of it form an orthonormal system in the Hilbert space 
of finite-energy signals. The latter pulse shape is described as the re- 
sponse of an infinite-order digital FIR filter on a Gaussian function 
as input signal. Several equivalent versions of the digital filter includ- 
ing their finite-order approximations are presented. All filters enjoy 
the property that explicit formulas for their coefficients and poles are 
availabele. The filters are fully parametrizable with respect to band- 
width and sampling rate of the digital data. 



1 Introduction 

The main new result of the paper is the signal described in Theorem [3] 
of Section 3. This signal is the root of an interpolating function found 
earlier in context with sampling theorems based on Gaussian functions [5], 
[6]. "Extracting the root" was prompted by [15J. The proof of Theorem 
[3] is given in some detail. It is based on spectral factorization and uses 
notions and identities from the theory of g-analogs [H] , [2] . Concepts of the 
latter theory normaly arise in the context of combinatorics, number theory 
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and statistical mechanics, in the field of communications they seem to be 
new (appart from coding theory where they have a different meaning). In 
Sections 4-5 new digital filters are presented. Background material on 
sampling in shift-invariant spaces is discussed in Section 2 in some length. 

The paper draws on previous work of Bhandari [3J, [3J, Unser [TU], [llj . 
ID, 03], Walter [8], and our own [5], [6], 0. 

The following notation is used: L 2 (R) is the space of square integrable 
functions (or finite-energy signals) / : R — > C with inner product (f\, /-j) = 
fi ( x )f2 {x) dx. For the Fourier transform we adopt the convention f(uj) = 
(27r) -1 / 2 e~ lxw fix) dx, where x denotes time and u (angular) frequency. 



2 Sampling in Shift-Invariant Spaces Revisited 

In the present section we give a review of sampling in shift-invariant spaces, 
its inherent difficulties and possible remedies. 

Suppose that ip E L 2 (IR) is a continuous function, satisfies ip(x) = 
0(\x\~ 1 ~ e ) as x — > ±oo for some e > and for any A E A, where A is some 
non-empty subset of M + = (0, +oo), the family of functions {</?(• — n\);n E 
Z} forms a Riesz basis in L 2 (R). Furthermore, we assume that for any A E A 

oo 

v(n\)e- inXuj + 0, cu E R. (1) 

n=— oo 

An example is the Gaussian function 

tp {x) = -=J e-^TW F rH T 0( w ) = -L e -£, ( 2 ) 

where the bandwidth parameter (5 is an arbitrary positive real number. 

The shift-invariant space V\{ip) is defined as the closed linear span of 
{</?(• — n\);n E Z} in L 2 (R). The following theorem asserts perfect re- 
construction of any / E V\(ip) from sample values /(nA), n E Z, for any 
sampling interval A E A. 

Theorem 1 Lei A E A. For any f E V\(v?) it holds that 

f(x) = J2f( nX ) c P^(x-n\),xeR, (3) 
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where the interpolating function ip- mt E V\(ip) is given 

£M yl = — 

'= EneZ + nA ) ' A 



This theorem is due to Walter [8], who proved it in case of an orthonormal 
system {</?(• — n);n G Z}. The argument carries over to Riesz bases {ip(- — 
nA);n G Z}, see [6] for a Gaussian function cp. Then, (fi nt is given by the 
right-hand side of Eq. (fTT|) after substitution of r by 2r (see below) [6j. 
Unser [13] uses shift-invariant spaces of the form 

„_ M ( X „\ „_ r- 



V T (<t>) = < f; f(x) = ^ or, n <t> [--n),ore£^,T>0, (4) 

as means for approximating functions / lying in other subspaces of L 2 (M), 
for instance in the Sobolev space W\. Then, for the vanishing of the ap- 
proximation error in the L 2 — norm as T — » it is necessary and sufficient 
that the generating function <p satisfies additionally the partition of unity 
condition (PUC), 

Vx G M : ^ <j)(x + n) = 1. (5) 

neZ 

In our approach, the shift-invariant space has the form [6] 



I n&L J 

{/,/(,) = EX= 



n I , c G £ 2 



where ^aO^) = Xip(Xx), i.e., now <j> = <p\. Since by the Poisson summation 
formula it holds that 

£ ^(x + n) = v^r £ (^) e 2 — , (6) 

neZ neZ ^ ' 

we retain an "approximate PUC" when A > becomes small enough under 
the assumption that |<^(w)| decays rapidly to as to — » ±oo and L ^(x) <ix = 



x We make no difference between ipi nt and y>i n t and so on. 
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1. Indeed, in case of a Gaussian function ([2]) we obtain with high accuracy 



^2ipx(x + n) = 1+ e"wF e 2 ™ K l,ieI, (7) 

neZ neZ\{0} 

as <- A < 2/p. 

Anyway, a serious problem arises when resampling elements of a shift- 
invariant space. In practice, one often has samples {/(a + nA)} of / G Vx(<p) 
with an offset a / 0. The following reconstruction formula can be derived 
(see [S] in case A = 1): 

/0) = J2^ a + nX )^t,ai x ~ nA )> » G M, (8) 

ragZ 

where the interpolating function <^i n t,a £ V\(p) is given by 



E„ez ¥>( a + nA)e" 



= ^ — | rv^Aj ' ( 9 ) 



(provided that the denominator does not vanish). Since <fint,a now depends 
on the often unknown offset a this raises a severe problem in practical ap- 
plications. 

One way out is to introduce a second function space, the so-called lo- 
calization space to be defined next. For the rest of the paper we restrict 
ourselves on a Gaussian function tp as given in ([2]), for the general approach 
see [TJ. Similarly to the classical sampling theorem we apply prior to sam- 
pling a prefilter 



(P 9 f)(x)= f{y)v(y-x)dy= e™"<p(u)f(u)ckj, (10) 

J — oo J — oo 

to an arbitrary finite-energy signal / G L 2 (1R) and the goal is to recover 
the filter output signal g = P v f from sample values g(nX), n G Z. The 
localization space of the prefilter, 

^ = { 5 = iV;/GL 2 (M)}, 

corresponds to the space of bandlimited signals in case of the classical sam- 
pling theorem. For the statement of the "approximate sampling theorem" 
(Theorem [2j see below) we also need the autocorrelation function of ip, 

$ = P v <p€ V v F r±? S( u ) = V^\0(u)\\ (11) 
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again a Gaussian function. For arbitrary A G A we define the (second) 
interpolating function <^i nt G 7^ by its Fourier transform 

*tatH = a „ ^ ; 71 = 27r/A - (12) 

Note that because of the Riesz basis condition on tp the denominator will 
never vanish [7]. In general #i n t 7^ -P^mt- 

Remark 1 By representation Iil2\) it is readily seen that 



A x 

& int (u; + nA) = l \/uj G R, (13) 



2vr 

n=— 00 

which is equivalent to <Pi n t(n\) = 5 n Vn G Z, a necessary and sufficient 
condition for the set of shifted pulses {$i n t(x — nX); n G Z} to 6e ISI-free at 
points in time x n = nX, n G Z J^]/. 

In the Gaussian case we compute that 

* intM = TiT^wvj' (14) 



where 



CO 



U^) = -L= £ e-v(^) 2 (15) 

\J —IT 

v n=— 00 

is one of the four theta functions pQ with parameter r given by 

r =JL(A/3) 2 (t = </=!). (16) 
By inversion of the Fourier transform we obtain 

1?^ (0, — -) smh(z7rrx/A) 

where i?i(z,r) = 2 ££1 g^+^V 1 )" sin[(2n + 1)ttz], q = e i7TT , 3(t) > 0, 
is again a theta function |T]. When A < it holds with high accuracy 
that ^int(^) ~ Sq{x) where 

f-y / \ . sin(7rx/A) ,„ „. 

S (x)=it , \ 1 18 
smh(«7rrx/ A) 

Note that in our setting it is always a negative real number. The following 
"approximate sampling theorem" is an instance of the Generalized Sampling 
Theorem [7J Theorem 3]. 
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Theorem 2 Let g be a signal in the localization space Vw of Pm, i.e., g = 
Ptpf for some function f G L 2 (R). Let A > and let the interpolating 
function <L> in t be given by ( f 17\ ). Then the function defined by 

oo 

~g{x)= g(n\)$ m t(x-n\),xem, (19) 

n=— oo 

is an element of the shift-invariant subspace lZ\((p) C 7? which may be 
defined as the closed linear span of {<?>(• — n\);n G Z} in L 2 (R). Perfect 
reconstruction is obtained at the sampling instants, 

g(n\) = g(n\), n G Z, 

while for any other i£l the reconstruction error is estimated by 

\g(x) - g(x)\ 2 < e~&- / \f(y)\ 2 dy. (20) 

V27T JjR 



Consequently, the reconstruction error decays superexponentially to as <— 
A < 1/0. 

In the specfific case of a Gaussian function as considered here it is shown 
in [6] that the reconstruction error becomes small already when A = 1//3 
(provided that the energy L |/(x)| 2 dx of the input signal / is not too large). 

Since the localization space is invariant with respect to arbitrary 
translations g —* g a = g(a + •), a G R, the reconstruction formula (fT9l) ap- 
plied to sample values {g(a + nA); n G Z} now simply results in an approx- 
imation g a G 1Z\(ip) to g a £ Vp. So, we get rid of the problems concerning 
an unknown offset mentioned earlier. 

Remark 2 The Fourier transform of function Sq as in U8\) is given by 

- , s A sinh(^) 

5b M = — 7= ? \ n \ ? A > 0, 

V ; cosh cosh 



and forms a simple though very good approximation to &i n t(cj), uj G R, when 
A < 1//3. 

3 Spectral Factorization of the Interpolating Func- 
tion 

We start by compiling some prerequisits. 
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Definition 1 For any a£l and g E C, \q\ < 1, the q-Pochhammer symbol 
(a; q) n is defined by 

(a;q) n = (1 - a)(l - aq) ■ ■ ■ (1 - ag*" 1 ),** = 1, 2, . . . 
(a;<?)oo = lim(a;g) n 

n^oo 

(a; g) = 1. 

The following identities of Euler (see, e.g., pj) hold for any q £ C, l^l < 1: 

00 fn(n-l) 00 

1 + E7^^ = IK 1 +^)'^ C ' ( 21 ) 

„=1 W,9Jn n=0 

OO n CO 

1 + E(i r = n( 1 -^r 1 ,N<i. (22) 

n=l ? ' n n=0 

Only Eq. (|21|) will be used in the proof of the next theorem. The Jacobi 
triple product identity [2] 

oo oo 

q n2 z n = ~[[(l-q 2n )(l+q 2n - 1 z- 1 )(l + q 2n ~ 1 z) (23) 

n=— oo n=l 

will also be needed. In our case always 

q = e™\ (24) 
where the parameter r is as in (|16j) . The special g-Pochhammer symbols^) 



{q 2 ;q 2 )n = \[{i-q 2k ),Qv := (g 2 ;g 
it=i 

will occur frequently. 
Theorem 3 For the signal 



2\ 

oo 



1 OO / \ r 

^ ortho(x) = 7f=2=2=, E (a 2. q n 2^ ^ X ~ nA )> ( 25 ) 

V w '9 Joo n=0 w 



v v - - - r. /(, — 

is i/ie Gaussian function (0|) normalized to unit energy, J K |</>(x)| 2 cte = 1, it 
holds that 

S- mt (u) = \/2^ ortho (u;)| 2 , uj G M. (27) 



2 In the theory of q-analogs Q3] the "g" is a generic name standing for q, q 2 , —q etc. 
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Proof. The definition (|14p of the interpolating function ^} nt in Fourier do- 
main may be written as 



Then, Eq. (|2"T|) becomes 

|Vortho(Wj| = A- 



Since the theta function $3(2, r) in fjl5[> has the second representation [T] 



00 

,2 



i? 3 (z,r) = l + 2^<f cos(2vrnz) 

n=0 

00 

= Yl q n2 e 2ninz , zeC, 

n=—co 

we obtain by setting 

00 

71=1 

in virtue of identity (|23p the factorization 



03(o;, t) = Q^ /2 P(e 2 ™) • Qj /2 P(e 2 ™), 

Since the real- valued function i?3(-,t) is positively lower bounded on M so 
is the function u> — > \P(e 2wtuJ )\. As a consequence, the definition in Fourier 
domain 

¥>orthoM = A2 - — — , 

(-ir) 1 /4gJ/ 2 p( e 2«u,/A) 

will result in a function i^ortho G £ 2 (M) satisfying condition ([27)1 . 

Next, we need to invert the Fourier transform. Since u — ► 1/P(e 27ria; ) 
is a bounded 1-periodic function it is in L 2 ([0, 1)) and has a Fourier series 
expansion 1/P(e 2 ™) = E^L-oc One" 2 ™", (a„)„ G z e £ 2 (Z), where 

le2 ™^^™^. (28) 



S 



By inverse Fourier transform we now obtain 
1 f°° ■ 1 

VorthoO) = C—= / e — — 1 - 7 -r-0(u) du, C 



V^J-oz P( e 2™M)^ ' (_ ir) i/4gi/2 

n= _oo V 2^ ./-oc 
00 

c a n cp(x — nX) 



n=— 00 

00 



E a n<t>ix-n\). 



^0 n=— 00 

The computation of the coefficients (|28p is accomplished in the complex 
domain. 

Case n = —1, —2, . . . : After substitution of uj by 1 — uj we obtain 

■l 1 



p(g-2™) 

1 r z - n - 1 

- 2^^1=1^) 

where integration in the latter integral is performed counterclockwise around 
the unit circle. Since the integrand function is analytic within a neighbour- 
hood of the closed unit disc, we obtain a n = 0, n = —1, —2, .... 
Case n — 0,1, ... : We now have 



1 f z n ~ l 



dz 



2m 7J Z | =1 P(z) 

lim (p — dz, 

M->oo 2m J\ z \ =l Pm{z) 

V V ' 



where 



M-l 



P M (z)= H(l+q 2k+1 z- 1 ),M = l,2,..., 

k=0 

and the path of integration is same as before. The integrand function in the 
integral defining Ijw(n) has simple poles in 



_2m+l 



, m = 0,1,... ,M- 1, 
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within the unit circle. By means of the theorem of residues we obtain 

M-i z n-i 
I A/ (n)=^Res 2mp - M , 



m=0 

whence 



00 n-1 

J " = JV M(n)= £ ReS -P^ 

m=0 v ' 



We compute that 

Res- ^— -r = Res 



z n-l 



P(z) m zP(z) 



(9W)ooTlT=i(i-<r 2fe ) 

^ ^jn / ^ jm^m(m+l) ^2n^m 

tewu nr=i(i-^) ' 

After substitution of q by g 2 Eq. (|2ip becomes 

00 n(n-l) 00 

1 + Er^r^ = n( 1 + ^ 2n )- 
^(<z 2 ;<? 2 )n ^ 



Hence 



(-?)' 



00 



7 m(m— 1) 



E « - 1 / „2(n+l) ym 



(<? 2 ;<? 2 )oc ^ (<? 2 ;<? : 



00 



(-g)" 



J^(l_ g 2(n + l )(? 2 fc) 



which concludes the proof of Theorem [3l 

Remark 3 (cf. |3j, [4j) a) It is well-known that for an arbitrary function 
s G L 2 (R) t/ie system of functions {s(x — n\);n S Z} forms an orthonomal 
system in I? (M) i/ and only if 

00 

A ^2 \K U + nA )\ 2 = 1 a - e -> 
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where A is as in IW) . For the signal s(x) = Vortho(^) in Theorem^ this is 
fulfilled because of Eqs. U3\) and \21^ (hence the subscript notation). 

b) At the receiver's side the filter P Voitho obtained after substitution in 
U0\) the Gaussian function ip by </? rtho forms a matched filter allowing opti- 
mal receipt of the signals s(x — nX) = ^p OT tho( x ~ i n the presence of noise. 
Moreover, because {v?ortho( - — nX);n £ 1*} forms an orthonormal system ISI 
is cancelled in the output of the receiver filter Pu, OItho at the points in time 
x n = nX, n S Z. 

In view of Theorem [2] we regard A = 1/(3 as the " natural" step size with 
respect to bandwidth parameter f3. 



4 Digital Filters 

Substitution of q by q 2 in Eqs. (i2~Tj) and (|22l) and, subsequently, of z by qz~ x 
and —qz^ 1 , resp., results in the identities 

oo n 2 oo 

n=l W ' q )n n=0 

°° (—n\ n °° 1 
l + Y X %■ z~ n = TT f-, r. (30) 

For N = 1,2,... we define the functions 



7V-1 



n=0 

N 



H?(z) = 1 + 53- 



^ fa 2 ;? 2 )* 

Since by Eqs. ([29]), ([30]) it holds that 



lim H?(z) = lim H?(z) = lim Izl > 1, 

for sufficiently large N the functions -£f^ (z), k = 1, 2, 3, are approximations 
of each other. 
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\z) is the transfer function of an IIR filter of order N, (z) of 
a FIR filter of the same order, and H2 (z) represents a cascade of ./V IIR 
filters of order 1. All three filters are causal filters. Recall that parameter 
q defined by ([To]) . (f!H|) depends on parameters j3 controlling bandwidth of 
the filter and A defining sampling interval of the discrete data. /3, A may 
be chosen independently and arbitrarily within (0, +00). H^(z) is a finite 
order approximation to the digital time-domain filter used in (|25|) . 

Because B-splines [12] with increasing order converge to Gaussian func- 
tions it is worthwile to compare digital filters based on Gaussian functions 
with their B-spline based counterparts |10j . [11] . The digital filter with 
transfer function 

N-l 1 
II (l + q 2n+l z -l)n +q 2n+l z \ 

H${z^)H?{z), 

a combination of a causal and an anticausal recursive filter, corresponds to 
the so-called direct B-spline filter [lOj (3.15)] occuring also in [31 Section 
IV]. Note that in contrast to the numerical values in [111 Table I] the poles 
of (z) within the unit circle are described by a simple formula, 

z n = -q 2n+1 ,n = 0,1,..., N-l, (31) 

where q € (0, 1) is explicitly given by (fT6j) . (j23T) . There are no poles on the 
unit circle. 



5 Implementation 

The digital filters with transfer functions (z), k = 1, 2, 3, of the previous 
section should now be implemented and evaluated. Note that the parameter 
q < 1 as in ([2^]) becomes close to 1 as the time-frequency product A/3 in (fl6|) 
becomes small. This would result in digital filters of high order N. Then, 
the recursive IIR filter H^(z) might be the best choice, because the factor 
n 2 in the exponent of the coefficients ensures faster decay resulting in a much 
lower order TV" needed. 

However, for the computation of the sample values V 9 ortho(w-A), m £ Z, 
of the presumably useful pulse form (j25j) by 

00 ( q) n 
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only a few filter coefficients are involved because of the fast decay of the 
Gaussian function <p as in (j26j) . 

A second useful pulse form is the interpolating function & mt of (JT3J) be- 
cause of its ISI-freeness. Here the question arises how to realize it technically. 
The solution given in form of the ISI-free kernel in [31 (7)] would consist in 
sending the Gaussian function <I> as in (jlip through a digital filter with co- 
efficients (^(A;))^ 1 (case A = 1). A finite order approximation is the digital 
filter H± (z) of the previous section (up to a constant factor we didn't care 
so far). The implementation of the related direct B-spline filter described 
in (2.5)] should carry over to our case with minor changes (if any). An 
implementation of the latter filter would be interesting on its own, see [3, 
Section IV] and [10, Section III-C]. 
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